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First-principles  kinetic  theory  is  used  in  this  work  to  analyze  non-equilibrium  homoge¬ 
neous  condensation  of  argon  and  water.  The  present  model  uses  a  recombination-reaction 
energy-dependent  mechanism  of  the  DSMC  method  for  the  dimer  formation,  and  RRK 
model  for  the  evaporation.  Three-step  validation  of  the  model  is  conducted,  (i)  comparison 
of  clusterization  rates  in  an  equilibrium  heat  bath  with  theoretical  predictions,  (ii)  compar¬ 
ison  of  the  argon  dimer  fractions  in  an  orifice  expansion  with  semi-analytical  correlations, 
and  (iii)  comparison  of  water  cluster  size  distributions  with  experimental  measurements. 
Reasonable  agreement  was  observed  for  all  three  parts  of  the  validation. 

I.  Introduction 

Two  different  approaches  for  modeling  the  condensation  in  rapidly  expanding  plumes  have  been  reported 
in  the  literature.  The  first  approach,  known  as  the  classical  approach,  takes  its  starting  point  from  the 
classical  nucleation  theory  (CNT)  which  is  based  on  equilibrium  thermodynamics.1,2  The  second  one,  known 
as  the  kinetic  approach,  treats  nucleation  as  the  process  of  kinetic  chemical  aggregation.3 

The  classical  approach  considers  the  energy  of  cluster  formation  from  the  vapor  state.  Assuming  unimolec- 
ular  reactions  of  cluster  growth  and  decay,  CNT  calculates  the  corresponding  condensation  and  evaporation 
rates  using  the  Gibbs  distributions  and  the  principle  of  detailed  balance. 3-6  The  nucleation  rate  is  then 
calculated  assuming  a  steady  state  condition.'  Although  a  rapidly  expanding  supersonic  plume  is  quite  dif¬ 
ferent  from  the  isothermal,  ideal  gas  environment  assumed  by  CNT,  the  classical  predictions  were  found  to  be 
qualitatively  correct  for  the  modeling  of  cluster  formation  in  supersonic  jets. 8-12  There  are  many  examples, 
however,  when  CNT-based  results  cannot  be  fitted  to  experimental  data.  The  CNT-based  prediction  of  the 
cluster  size  distributions13  significantly  deviated  from  experimental  data.14  This  is  consistent  with  the  work 
of  Ref.  15,  where  it  was  also  found  that  the  correct  prediction  of  the  cluster  size  distribution  along  with  the 
internal  and  translational  energy  distributions  is  beyond  the  area  of  applicability  of  the  classical  approach. 

The  reasons  for  the  discrepancy  between  the  CNT-based  distributions  and  experimental  data  are  both  due 
to  problems  inherent  in  CNT  and  the  flow  conditions  of  expanding  plumes.  The  former  include  the  ambiguous 
definition  of  the  surface  energy  of  small  clusters,7  the  negligence  of  the  rotational  and  translational  degrees 
of  freedom  of  freshly  nucleated  clusters,16  and  the  unrealistic  description  of  vapor-cluster  and  cluster-cluster 
interactions.12  The  latter  are  related  to  the  main  assumptions  underlying  the  derivation  of  the  nucleation 
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rate,  which  may  be  violated  in  rapidly  expanding  supersonic  flows.17  The  transient  time  needed  for  a  system 
to  reach  steady  state  in  terms  of  the  unimolecular  cluster  reactions  may  be  such  that  the  jet  macroparameters 
will  significantly  change  during  that  time.  Moreover,  many  theoretical  and  experimental  results18-20  suggest 
that  local  thermal  equilibrium  does  not  exist  in  an  expanding  supersonic  jet.  Thus  the  process  of  cluster 
formation  is  not  likely  to  be  isothermal.  There  have  been  recent  advances  in  CNT  mainly  aimed  at  achieving 
a  more  realistic  model  for  condensation  and  evaporation  rates,6’ 16,21-24  but  other  principal  deficiencies  of 
CNT  and  its  application  to  the  non-equilibrium  environment  still  have  yet  to  be  addressed.  Note  that  despite 
these  deficiencies,  the  steady  state  CNT  nucleation  rate  with  a  correction6  is  used  in  such  commercial  fluid 
dynamics  codes  as  GASP  or  CFD-particle. 

Unlike  CNT,  the  kinetic  approach  does  not  assume  local  thermodynamic  equilibrium.  Instead,  a  micro¬ 
scopic  view  of  the  interactions  of  monomers  and  clusters  is  established  either  analytically  via  a  mathematical 
model,  e.g.,  by  the  Smolucliowski  equations  where  the  interaction  between  particles  is  modeled  by  the  reaction 
rates,25-27  or  in  computer  simulations,  e.g.,  in  molecular  dynamics  (MD)  calculations  where  the  interaction  is 
modeled  by  an  interaction  potential.28-30  It  can  be  shown  that  the  application  of  the  Smoluchowski  equations 
to  the  modeling  of  clustering  in  supersonic  jets  is  computationally  unfeasible.  Even  though  MD  simulations 
might  seem  attractive  (since  no  information  besides  the  interaction  potential  is  needed  to  perform  the  calcu¬ 
lations)  ,  they  are  computationally  limited  to  a  system  size  of  about  one  million  gas  particles  and  a  time  scale 
of  a  few  nanoseconds  at  most.  In  real  plumes,  either  thruster  nozzle  or  ablation  originated,  the  number  of  gas 
particles  and  the  expansion  time  is  greater  by  many  orders  of  magnitude  (see,  for  example,9,31).  Therefore, 
the  MD  technique  cannot  be  directly  applied  to  the  simulation  of  even  small  laboratory-sized  supersonic 
jets. 

A  promising  direction  in  modeling  the  coupled  condensation  flow  is  the  use  of  a  kinetic  particle  simulation 
method,  direct  simulation  Monte  Carlo  (DSMC),32  which  is  applicable  in  a  wide  range  of  flow  regimes  from 
free  molecular  to  near  continuum.  It  is  a  statistical  approach  for  solving  the  spatially  nonuniform  master 
Leontovich  equation  for  the  N-particle  distribution  function  and,  in  the  limit  of  a  large  N,  it  represents  an 
accurate  solution  of  the  Boltzmann  equation.  The  advantage  of  DSMC  as  compared  to  other  methods  is  that 
complicated  cluster-cluster  and  cluster-monomer  interactions  including  the  multi-body  reactions  of  cluster 
nucleation  can  be  seamlessly  incorporated. 

The  DSMC  method  has  been  used  to  study  the  process  of  cluster  formation  and  evolution  for  a  number 
of  years.5,33,34  However,  the  gas  flow  in  the  earlier  studies  was  uniform,  the  considered  cluster  size  range  was 
very  narrow  (up  to  25  monomers  in  a  cluster)  and  the  examined  reaction  types  were  unrealistically  limited 
to  elastic  collisions,  cluster  and  monomer  sticking  to  clusters,  and  evaporation  of  monomers  from  clusters. 

More  recently,  the  DSMC  method  has  been  extensively  and  successfully  applied  to  modeling  the  processes 
of  cluster  formation  and  evolution  in  supersonic  jets  by  Levin  et  al  (see,  for  example,  Refs.  12,35,36).  The 
model  initially  was  based  on  the  classical  nucleation  theory,  with  the  new  clusters  being  formed  at  the  critical 
size.  Further  work  of  these  authors37  extended  the  kinetic  dimer  formation  approach  of  Ref.  38,  who  assumed 
that  a  ternary  collision  always  results  in  a  dimer  formation,  to  include  molecular  dynamic  (MD)  simulations 
for  obtaining  information  on  the  probability  of  dimer  formation  in  such  ternary  collisions  The  work39  used 
a  temperature-dependent  probability  of  formation  of  argon  dimers. 

Although  the  use  of  MD  has  a  number  of  advantages,  due  to  its  inherent  limitations,  the  obtained 
probability  cannot  be  unambiguously  related  to  such  characteristics  of  the  ternary  collision  as  the  internal 
energy  of  the  collision  complex  and  the  kinetic  energy  of  the  impinging  monomer.  Another  possible  limitation 
of  the  above  work  is  that  even  though  the  nucleation  process  is  fully  kinetic  in  the  most  recent  papers,  they 
the  CNT  rate  is  used  for  evaporation,  although  there  may  be  more  sophisticated  models  for  nonequilibrium 
cluster  physics.40 

In  the  present  paper,  the  DSMC  approach  for  modeling  of  homogeneous  nucleation  in  rapidly  expanding 
plumes  is  extended  to  include  a  number  of  new  features.  Most  importantly,  a  truly  kinetic  RRK  model41  is 
implemented  to  characterize  the  cluster  evaporation  rates.  Then,  an  energy  dependent  collision  procedure 
similar  to  the  recombination  reaction  model  of  Ref.  42  is  used  for  the  collision  complex  formation.  An 
empirical  parameter  is  used  for  the  inelastic  collision  number  in  the  cluster-monomer  collisions.  For  dimers, 
this  parameter  was  calibrated  through  the  comparison  of  the  computed  nucleation  rates  and  equilibrium 
constants  in  thermal  bath  with  available  theoretical  and  experimental  data  for  argon  and  water.  Addi¬ 
tional  validation  analysis  is  conducted  through  comparison  with  the  empirical  correlation43  for  argon  and 
measurements  of  cluster  size  distribution  in  water  expansions.44 
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II.  Computational  model  of  homogeneous  nucleation 


In  the  present  first-principle  model  of  homogeneous  condensation  formulated  for  the  DSMC  method, 
all  of  the  most  important  processes  of  cluster  nucleation  and  evolution  are  considered  at  the  microscopic 
level.  The  main  processes  that  are  included  in  the  model  and  described  in  detail  below  are  (i)  formation  of 
collision  complexes  through  the  binary  collisions  of  cluster-forming  monomer  species,  (ii)  creation  of  dimers 
through  the  collision  stabilization  of  collision  complexes,  (iii)  elastic  monomer-cluster  collisions  that  change 
the  translational  and  internal  energies  of  colliding  particles,  (iv)  inelastic  monomer-cluster  collisions  that 
result  in  monomer  sticking,  (v)  cluster-cluster  coalescence,  (vi)  evaporation  of  monomers  from  clusters.  The 
details  on  each  of  these  processes  are  given  below. 


A.  Collision  complex  formation  and  stabilization 

One  of  the  important  assumptions  of  the  present  model  is  that  all  pairs  of  colliding  particles  create  collision 
complexes.  A  collision  complex  is  a  pair  of  monomers  that  have  collided,  and  may  have  the  conditions 
necessary  to  form  a  dimer  if  struck  by  a  third  particle  during  its  lifetime.  The  collision  complex  lifetime,  t; , 
is  assumed  to  be  dependent  on  the  type  of  monomers  and  their  relative  collision  velocity,  with  the  functional 
dependence  given  by  the  well  known  Bunker’s  expression45 

ti  =  1.5ao^  eft  E~i,  (1) 


where  a o  and  eo  are  the  potential  depth  and  separation  distance  parameters  of  the  Lennard- Jones  potential, 
g  is  the  reduced  mass  of  the  colliding  particles,  and  E  is  their  relative  translational  energy.  The  values  of 
do  and  e0  used  in  this  work  are  3.166  x  10-10  m  and  1.079  x  10-21  J  for  water,  and  3.405  x  10_1°  m  and 
1.654  x  10-21  J  for  argon. 

The  process  of  interaction  of  collision  complexes  with  surrounding  gas  particles  is  modeled  using  the 
majorant  frequency  scheme46  with  the  assumption  that  the  collision  complex  -  third  particle  interactions 
are  governed  by  the  Variable  Hard  Sphere  (VHS)  interaction  model.47  The  VHS  parameters  of  a  collision 
complex  were  assumed  to  be  those  of  the  comprising  monomers. 

Generally,  the  probability  P(t)  that  a  collision  complex  will  collide  with  a  third  particle  during  an 
arbitrary  time  t  is 


dP(r) 

dr 


=  v{1-P{t)) 


(2) 


where  v  is  the  collision  frequency  of  the  collision  complex  with  third  particles.  For  a  mixture  of  Ns  gas 
species,  v  is  expressed  as 

Na 

v  =  'Y^ni((jg),  (3) 


where  a  is  the  corresponding  total  collision  cross-section,  g  is  the  relative  collision  velocity  between  particles, 
and  brackets  denote  averaging  over  g.  Obviously,  the  probability  that  the  collision  complex  will  move  freely 
and  not  collide  during  r  is  given  by 


Pfree{^~ )  —  & 


—  A— jG 


(4) 


This  expression  represents  the  probability  that  no  dimer  will  be  formed  during  r.  Using  the  inverse  transform, 


tc  =  -  In {R/v), 


(5) 


where  R  is  a  random  number  uniformly  distributed  between  0  and  1.  The  majorant  frequency  algorithm  for 
each  pair  of  colliding  monomers  is  therefore  as  follows. 

1.  Calculate  t/  and  Vmajorant  and  set  t0  =  0,  *  =  1.  Here,  vmajorant  =  namaXgmax  is  the  majorant  collision 
frequency,  and  amax  and  gmax  are  maximum  cross  section  and  relative  collision  velocity,  respectively. 

2.  Calculate 

U=ti- i-ln( - — - j  (6) 

\  ^ majorant  J 

3.  If  ti  >  At,  where  At  is  the  simulation  timestep,  go  to  the  the  next  pair  of  monomers. 

4.  If  ti  <  At,  then  a  physical  collision  occurs  with  a  probability  Pc  =  JmaxgmaX ■  With  a  probability 
1  —  Pc,  i  =  i  +  1  and  the  algorithm  returns  to  step  2. 
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In  this  algorithm,  for  consistency  with  the  collision  complex  lifetime  determination,  an  expression  for  the 
diameter  d  of  the  collision  complex  recommended  in  Ref.  45  was  used, 

d  =  31ia0(e0/E)i.  (7) 


B.  Dimer  stabilization 


If  there  is  a  physical  collision  between  a  collision  complex  and  a  third  particle,  there  is  a  possibility  of  forming 
a  stable  dimer  as  a  result  of  such  a  collision.  Generally,  the  probability  of  the  formation  of  a  stable  dimer 
(dimer  stabilization)  depends  on  the  colliding  species  and  the  energies  -  both  translational  and  internal  -  of 
the  colliding  particles.  In  this  work,  constant  stabilization  probabilities  of  0.2  for  Ar37  and  1.0  for  H2O  were 
assumed,  which  seem  reasonable  for  the  range  of  temperatures  under  consideration. 

Dimer  creation  through  the  collisional  stabilization  of  collision  complexes  is  a  two-step  process,  L  +  M  —* 
(LAI),  (LAI)  +  K  — *  LAI  +  K.  Here,  L  and  AI  are  monomers,  (LAI)  is  the  collision  complex,  and  K  is  the 
third  particle.  The  algorithm  of  this  process  is  described  below. 

1.  Velocities  of  the  collision  complex  are  calculated  from  the  momentum  conservation  as 


V(LM)  = 


mLvL  +  mMVM 

^(IM) 


(8) 


2.  The  internal  energy  of  the  collision  complex  is  calculated  from  energy  conservation, 


EitfM)=EiLnt  +  Et}t  +  AE,  (9) 

where  A E  =  m.(LM)V2LM)  -  mKv2L  -  mLv2M. 

3.  The  total  energy  of  the  collision  complex  -  third  particle  pair  is  increased  by  the  evaporation  (or  binding) 
energy  Eevap, 

rptotal  rprel  i  rpint  ,  Tpint  ,  rp 

&  J-J(LM)  —  K  '  ^(LM)  '  +  &evap- 

Here,  E^M\_K  is  the  relative  translation  energy  of  the  (LAI)  —  K  collision.  Evaporation  energy  is  a  function 
of  cluster  size,  and  the  values  used  for  Ar  and  H2O  are  given  in  the  following  sections.  Note  that  E 'J?4  may 
be  used  instead  of  E^1,  since  the  vibrational  mode  of  third  particles  will  barely  be  excited  at  the  low 
temperatures  at  which  homogeneous  condensation  usually  occurs. 

4.  New  energies  E^lM)_K,  E are  sampled  using  the  Larsen-Borgnakke  scheme48  extended  to 

multiple  energy  modes.32 


C.  Elastic  and  inelastic  reflective  collisions  of  monomers  and  clusters 

The  interactions  between  monomers  and  clusters  is  one  of  the  key  processes  that  determine  the  nucleation 
rate.  The  reason  for  it  is  strong  dependence  of  the  evaporation  rate  on  the  cluster  internal  energy.  Since 
the  monomers  are  dominant  in  the  flows  considered  here,  the  cluster  internal  energy  is  mostly  governed  by 
its  relaxation  through  cluster-monomer  collisions.  In  this  work,  a  hard  sphere  model  is  assumed  for  cluster- 
monomer  collisions,  with  the  cluster  diameter  determined  from  Eqn.  7  for  dimers,  and  for  larger  clusters 
from  an  empirical  correlation  used  extensively  in  the  past  (see,  for  example,  Ref.  12), 

d=  2-  (A-ii  +B),  (11) 

where  A  and  B  are  species-dependent  constants,  and  i  is  the  number  of  monomers  in  the  cluster.  In  this 
work,  the  values  of  A  and  B  were  2.3  x  10~10  m  and  3.4  x  10”10  m  for  argon,  and  1.9  x  10_1°  m  and 
2.4  x  10_1°  m  for  water. 

For  the  energy  transfer  between  the  relative  translational  and  internal  modes  of  the  cluster  and  monomer, 
the  Larsen-Borgnakke  model  is  used,  and  a  parameter  Z  is  introduced  that  has  a  meaning  of  the  internal 
energy  relaxation  number.  The  energy  transfer  between  all  energy  modes  of  the  cluster-monomer  pair  occurs 
with  a  probability  Z  ,  and  an  elastic  collision  with  no  internal  energy  exchange  occurs  with  the  additional 
probability  1  —  Z~x .  A  value  of  Z  =  6  was  used  in  the  present  calculations,  that  allows  reasonable  agreement 
of  the  DSMC  rates  for  dimer  nucleation  and  dissociation  with  rates  available  in  the  literature. 
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D.  Sticking  collisions  of  monomers  and  clusters 


When  a  monomer  collides  with  a  cluster,  sticking  of  the  monomer  to  a  cluster  surface  is  possible,  in  addition 
to  a  reflective  collision  described  in  the  previous  section.  For  small  clusters,  monomer  sticking  is  the  main 
process  that  governs  the  evolution  of  the  droplet  size  distribution.44  An  empirical  dependence  of  the  sticking 
probability  on  the  species  radius  and  mass,  given  in  Ref.  49,  is  used  in  this  work.  This  dependence  reduces 
to 

e  =  dl  f  mn 

d 2  +  d\  \mn  +  mi 


where  indices  n  and  1  refer  to  the  cluster  and  monomer,  respectively. 

The  following  algorithm  is  used  to  model  sticking  of  a  molecule  to  a  cluster  L  to  form  a  larger  cluster  K. 
1.  Calculate  velocities  of  K  as 


v  k 


mLVL  +  mMVM 
mk 


(13) 


2.  After  velocities  vk  are  assigned,  the  new  internal  energy  is  calculated 


Eg*  =  Eg*  +  Eg*  -Q  +  AE 


(14) 


Here, 

A  E  =  -mKv%  +  mLv2L  +  mMv2M 

and  Q  =  Eevap  is  the  evaporation  energy  of  one  monomer  off  cluster  K. 

The  only  outcome  for  cluster-cluster  collisions  is  assumed  to  be  coalescence.  The  algorithm  for  cluster- 
cluster  sticking  collisions  is  similar  to  the  monomer-cluster  collision,  with  the  exception  of  Q  =  — Qk  +  Ql  + 
Qm,  where  Qi  is  the  energy  of  vaporization  of  cluster  i. 


E.  Evaporation  rate  and  algorithm 


Following  Ref.  31,  RRK  theory  is  used  to  model  the  evaporation  process.  The  evaporation  rate  ke  is  calculated 
as 


ke 


vNs 


-  E 
Eint 


evap 


3n—7 


(15) 


Here,  n  is  the  number  of  monomers  in  the  cluster,  v  is  the  vibration  frequency,  N s  is  the  number  of  surface 
atoms,  and  Eint  is  the  cluster  internal  energy.  For  dimers,  the  the  exponent  3n  —  7  is  replaced  with  1.  The 
number  of  surface  atoms  is  n  for  N  <  5,  n  —  1  for  4  <  n  <  7,  and  (307r)1/3(n1/3  —  l)2  for  n  >  6.  The 
vibration  frequency  was  taken  to  be  2.68  x  1012  s_1  for  water  clusters,  and  1012  s_1  for  argon  clusters.  Since 


dNs 

dt 


- keNs , 


(16) 


the  time  r  to  the  next  evaporation  event  for  a  given  cluster  may  be  sampled  from  r  =  —  In  .  The  algorithm 
used  to  model  the  evaporation  process  for  a  given  cluster  over  a  simulation  timestep  AT  is  as  follows. 

1.  Set  tlocal  d 

2.  Calculate  ke 

3.  Change  tiocai  —  tlocal— l  /n(R) / ke 

4.  If  tiOCai  >  At,  exit. 

5.  Evaporate  one  monomer  (see  below). 

6.  If  the  remaining  cluster  is  a  monomer,  exit.  Else,  go  to  Step  3. 


F.  Energy  redistribution  in  evaporation 

The  following  energy  redistribution  scheme  is  used  to  model  the  evaporation  of  a  monomer  M  off  a  cluster 
K,  with  a  smaller  cluster  L  formed. 

1.  Decrease  the  internal  energy  of  K  by  the  evaporation  energy  Eevap 

rpint  _  rpint  rp  /i7\ 

■t^K  ~  ~  ^evap  {*- •  ) 


5  of  16 

American  Institute  of  Aeronautics  and  Astronautics  Paper  2009-3745 
Distribution  A:  Approved  for  public  release;  distribution  unlimited 


2.  Split  cluster  internal  energy  E™*  between  the  relative  translational  energy  Erel  and  the  internal  energy 
of  the  M  L  pair  (the  sum  of  the  smaller  cluster  internal  energy  and  the  monomer  internal  energy),  Eint 
using  the  Larsen-Borgnakke  procedure. 

3.  Split  energy  Eint  between  the  smaller  cluster  internal  energy  and  the  monomer  internal  energy  using  the 
Larsen-Borgnakke  procedure. 

4.  Calculate  new  smaller  cluster  and  monomer  velocities  using  Erei ,  keeping  in  mind  that  the  center  of  mass 
of  the  new  pair  has  the  same  velocity  as  the  velocity  of  the  original  cluster. 

G.  Evaporation  energy  and  the  number  of  internal  degrees  of  freedom 

It  is  important  to  use  reasonable  values  for  the  evaporation  energy  of  a  monomer  off  a  cluster  and  the 
number  of  cluster  internal  degrees  of  freedom  as  a  function  of  cluster  size.  The  values  of  the  evaporation 
energy  where  taken  from  Refs.  50,51  for  water  clusters  and  Ref.  52  for  argon  clusters.  The  number  of  cluster 
internal  degrees  of  freedom  is  calculated  from  the  expression  for  the  average  internal  energy  (E)  of  a  cluster 
of  a  size  n 

eint  o 

(E)  =  ^kT  =  nCvT--kT 
as 

ent  =  n2^  _  3 
k 

where  Cv  is  the  cluster  heat  capacity.  For  water,  the  values  of  Cv  where  taken  from  Refs.  50,51,  while  for 
argon,  they  where  assumed  to  approximate  an  expression  =  2(3n  —  77),  where  77  is  5  for  n  =  2  and  6 
otherwise.53  Both  the  heat  capacity  and  evaporation  energies  are  listed  in  Table  G.  For  the  cluster  sizes 
larger  than  given  below,  the  values  for  the  maximum  listed  sizes  are  used. 


Cluster  size 

H2O  Eevap,  J 

h2oc„,  j/k 

Ar  Eevap,  J 

Ar  Cv,  J/K 

2 

4.919E-20 

8.347E-23 

0.989E-21 

0.17E-22 

3 

6.329E-20 

9.773E-23 

1.978E-21 

0.21E-22 

4 

1.021E-19 

1.029E-22 

2.967E-21 

0.26E-22 

5 

1.122E-19 

1.070E-22 

3.84E-21 

0.29E-22 

6 

1.009E-19 

1.106E-22 

4.48E-21 

0.31E-22 

7 

7.210E-20 

1.138E-22 

5.92E-21 

0.33E-22 

8 

6.889E-20 

1.166E-22 

4.70E-21 

0.34E-22 

9 

6.729E-20 

1.192E-22 

6.08E-21 

0.34E-22 

10 

6.569E-20 

1.218E-22 

6.13E-21 

0.35E-22 

11 

6.409E-20 

1.234E-22 

6.24E-21 

0.36E-22 

12 

6.409E-20 

1.243E-22 

6.64E-21 

0.36E-22 

13 

6.409E-20 

1.250E-22 

5.76E-21 

0.37E-22 

Table  1.  The  water  and  argon  cluster  heat  capacities  and  evaporation  energies  per  monomer. 


The  first  principle  condensation  model  described  here  was  implemented  in  the  DSMC  code  SMILE.54 
The  validation  of  the  code  through  the  comparison  with  theoretical  and  experimental  results  is  presented 
below. 


III.  Thermal  bath  relaxation 

Inelastic  cross  sections  for  monomer-monomer  and  monomer-cluster  collision  processes  are  necessary  for 
detailed  validation  of  a  kinetic  condensation  model.  These  cross  sections,  generally  functions  of  the  energy 
states,  both  translational  and  internal,  of  pre-  and  post-collisional  particles,  are  not  available  for  most  gas  and 
temperature  conditions  of  interest.  Contrary  to  the  energy  dependent  cross  section,  the  integral  temperature 
dependent  rates  for  such  collisions  at  conditions  close  to  equilibrium  are  available  in  the  literature.  Therefore, 
one  of  the  key  indicators  of  the  accuracy  and  reliability  of  a  condensation  model  is  its  ability  to  produce 
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realistic  rates  of  evaporation  and  nucleation  at  equilibrium.  Although  matching  the  rates  generally  does  not 
guarantee  correct  behavior  in  nonequilibrium,  it  still  is  a  necessary  condition  for  a  model  to  satisfy. 

In  this  work,  thermal  bath  relaxation  of  pure  argon  and  pure  water  are  examined  at  different  temperature 
conditions,  and  the  dimer  formation  rates  for  argon  and  equilibrium  constants  for  the  formation  of  dimers  in 
argon  and  water  are  calculated  and  compared  to  the  published  results. 55-57  In  all  thermal  bath  results,  one 
million  simulated  particles  were  used,  and  the  run  proceeded  until  the  steady  state  is  reached,  after  which 
the  results  were  sampled  for  20  thousand  timesteps.  The  number  density  was  5  x  1023  molec/m3  for  argon 
and  2  x  1023  molec/m3  for  water.  The  timestep  of  2.5  x  10-11  s  was  selected  so  that  the  number  of  collisions 
per  molecule  is  much  smaller  than  unity,  and  the  results  are  independent  on  the  timestep. 

The  computed  dimer  formation  rates  krec  for  argon  are  presented  in  Fig.  1,  where  they  are  compared 
with  the  stable  dimer  formation  rate  of  Ref.  55,  where  they  were  calculated  using  classical  trajectories,  and 
the  following  expression  was  proposed, 

krec  =  10.15 T~0-278  exp  {— 0.0031T}  . 

Generally,  the  present  dimer  formation  rates  are  in  reasonable  agreement  with  the  classical  trajectory  cal¬ 
culations,  with  the  maximum  difference  approaching  20%  for  higher  temperatures.  Note  that  the  computed 
rate  has  visibly  smaller  slope  than  that  of  Eqn.  III.  A  number  of  factors  could  be  affecting  the  slope,  among 
them  are  the  temperature  dependence  of  the  dimer  stabilization  probability  and  the  heat  capacity,  that  were 
not  included  in  the  present  model,  as  well  as  the  evaporation  energy  value  and  the  after-reaction  energy 
redistribution  (recall  that  the  Larsen-Borgnakke  model  was  used  for  the  energy  redistribution). 


Figure  1. 


Argon  dimer  formation  rate  and  equilibrium  constant  as  functions  of  gas  temperature. 


The  computed  equilibrium  constant,  Keq  that  is  the  ratio  of  the  dimer  dissociation  to  the  dimer  formation 
rate,  is  also  given  in  Fig.  1.  It  is  compared  to  the  theoretical  results  of  Ref.  56,  where  a  number  of  approximate 
classical  and  quantum  method  are  compared  with  exact  numerical  calculations.  Note  that  while  the  results 
for  different  models  and  interaction  potentials  where  found  to  be  widely  different,  there  was  a  good  agreement 
between  analogous  quantum  and  classical  calculations.  Figure  1  shows  that  there  is  a  very  good  agreement 
for  temperatures  of  200  K  and  higher,  whereas  for  100  K  the  present  rate  underpredicts  that  of  Ref.  56. 
Similar  to  the  dimer  formation  rates,  the  slope  for  the  equilibrium  constant  obtained  with  the  present  model 
is  smaller  than  the  theoretical  one;  the  reasons  for  this  may  be  similar  to  those  listed  for  the  dimer  formation 
rates. 
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An  important  factor  that  influences  the  magnitude  of  the  equilibrium  constant  Keq  was  found  to  be  the 
probability  of  the  energy  transfer  between  the  internal  modes  of  a  dimer  and  the  translational  modes  in 
dimer- monomer  collisions,  .  This  may  be  explained  as  follows.  The  dimers  are  formed  after  three-body 

collisions,  and  typically  have  internal  energies  smaller  than  the  evaporation  energy  after  those  collisions.  In 
argon,  the  evaporation  energy  for  a  dimer  is  relatively  small  compared  to  the  typical  total  collision  energy 
for  all  temperatures  under  consideration  ( Eevap/k  ss  70  K).  That  means  that  most  of  the  dimers  will  have 
their  internal  energy  in  excess  of  the  evaporation  energy  just  after  one  or  two  collisions  with  monomers. 
The  lifetime  of  the  dimers  whose  internal  energy  is  larger  than  the  evaporation  energy  is  very  short,  on  the 
order  of  a  picosecond.  This  results  in  the  dimer-monomer  energy  transfer  being  the  main  process  that  leads 
to  quick  dimer  dissociation.  Note  that  the  value  of  Z  has  negligible  impact  on  the  dimer  formation  rates, 
and  only  the  evaporation  rates  are  affected.  As  a  result,  in  the  range  of  temperatures  considered  in  this 
work,  the  equilibrium  constant  for  argon  was  found  to  be  nearly  proportional  to  Z~x .  For  example,  at  200  K 
Keq  =  0.11  x  10  26  m3/molec  when  Z  =  1  is  used,  versus  Keq  =  0.66  x  10  26  m3/nrolec  for  the  baseline 
case  of  Z  =  6. 

The  Z  dependence  of  the  equilibrium  constant  is  quite  different  for  water  molecule  condensation.  In  this 
case,  the  evaporation  energy  of  a  dimer  is  much  larger  than  the  translational  energy  of  colliding  molecules 
and  dimers  (the  reduced  evaporation  energy  Eevap/k  «  3, 500  K,  compared  to  gas  temperatures  on  the  order 
of  300  K).  The  high  value  of  the  evaporation  energy  results  in  longer  lifetimes  of  dimers,  since  much  more 
collisions  are  necessary  to  transfer  enough  energy  from  the  translational  modes  to  the  internal  modes  of  a 
dimer.  The  dependence  of  Keq  on  Z  is  therefore  much  weaker  for  water  than  for  argon.  In  a  250  K  thermal 
bath,  Keq  decreases  from  0.605~26  m3/molec  to  0.343  x  10~26  m3/molec  when  Z  decreases  from  6  to  1.  Note 
that  since  the  number  of  internal  degrees  of  a  dimer  is  large,  typically  about  22,  the  average  reduced  energy 
per  mode  is  ~160  K.  This  results  in  the  energy  in  dimer-monomer  collisions  being  mostly  transferred  from 
the  translational  modes  to  the  dimer  internal  modes  for  the  gas  temperatures  higher  than  160  K  (which  is 
the  case  in  most  of  the  flow  scenarios  considered  in  this  work). 

Comparison  of  the  equilibrium  constant  obtained  with  the  present  model,  with  the  theoretical  results  of 
Ref.  57,  where  a  flexible  potential  energy  surface  fitted  to  spectroscopical  data  was  used,  is  shown  in  Fig.  2. 
Note  that  there  are  a  number  of  theoretical  predictions  of  the  equilibrium  rate  of  water  dimerization,  and 
they  differ  by  at  least  a  factor  of  three  in  the  range  of  temperatures  considered  in  this  work.  Reference  57  was 
chosen  for  comparison  as  the  most  sophisticated  and  one  of  the  most  recent  ones.  The  following  expression, 
adapted  from  Ref.  57,  is  used  here, 

Keq  =  6.5202  x  10~32Texp  {1851.09/T  -  5.10485  x  1(T3T}  . 

There  is  a  very  good  agreement  for  gas  temperatures  above  270  K.  For  lower  temperatures,  the  results 
obtained  with  the  present  model  start  to  deviate.  The  difference  reaches  60%  at  200  K,  with  the  present 
results  being  smaller.  This  general  trend  of  underprediction  at  lower  temperatures  is  similar  to  that  of  argon 
bath. 


IV.  Dimer  mole  fraction  correlation  in  argon  orifice  expansion 

An  additional  key  measure  of  the  accuracy  and  reliability  of  a  condensation  model  is  the  terminal  dimer 
mole  fraction  from  a  sonic  orifice,  where  a  semi-empirical  correlation  is  available,43  based  on  a  large  number  of 
experimental  data  points.  In  this  work,  argon  expansion  from  a  sonic  orifice  is  modeled  at  several  different 
source  temperatures  and  number  densities.  The  computational  domain  consists  of  a  section  of  a  plenum 
2.5  mm  in  radius  and  2.5  mm  in  length,  along  with  an  expansion  region  that  is  2.5  mm  in  radius  and  7  mm 
in  length.  The  DSMC  runs  use  on  average  about  10  million  simulated  molecules  and  2  million  collision  cells. 
The  terminal  dimer  mole  fractions  are  recorded,  averaged,  and  compared  to  published  results.43 

Generally,  the  terminal  dimer  mole  fraction  occurs  when  the  expanding  gas  is  at  low  enough  densities 
to  be  nearly  collisionless,  and  the  cluster  mole  fraction  reaches  an  equilibrium  value.  As  seen  in  Fig.  3, 
this  occurs  well  within  the  computational  domain.  4  mm  downstream  of  the  orifice,  number  density  in  a 
typical  case  drops  by  2  orders  of  magnitude.  Near  the  outflow  boundary,  the  number  density  decreases  by 
yet  another  order  of  magnitude,  and  the  gas  mean  free  path  in  that  region  is  on  the  order  of  the  orifice 
diameter.  Note  that  the  figure  is  stretched  in  the  radial  direction  to  allow  better  detail. 

Cluster  mole  fraction  reached  a  terminal  value  in  all  of  the  cases  considered  in  this  work,  as  seen  in  Fig.  4 
(left).  Note  that  both  higher  number  densities  and  lower  temperatures  increase  both  source  and  terminal 
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Temperature  (K) 

Figure  2.  Water  equilibrium  constant  as  a  function  of  gas  temperature. 
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Figure  3.  Number  density  flowfield  from  a  sonic  orifice  with  a  source  temperature  and  number  density  of 
200  K  and  1.2  x!024  m-3,  respectively. 
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cluster  mole  fractions.  It  is  also  worth  noting  that  the  mole  fractions  in  the  plenum  are  in  good  agreement 
with  the  equilibrium  constants  calculated  in  Ref.  56.  This  figure  is  also  indicative  of  mole  fraction  of  dimers 
along  the  plume,  as  the  clusters  are  mainly  dimers  with  a  small  number  of  trimers  at  most  points  in  the 
plume. 


Figure  4.  Cluster  mole  fractions  along  the  axis  for  all  cases  considered,  left,  and  average  terminal  dimer  mole 
fraction  along  the  axis  as  a  function  of  scaling  parameter  SP. 

Finally,  the  average  terminal  dimer  mole  fraction  is  shown  in  Fig.  4  (right)  as  a  function  of  a  scaling 
parameter  SP.  SP  is  given  in  Ref.  43  as:  SP  =  n0 The  results  of  this  work  are  shown 

compared  to  an  empirical  expression  derived  in  Ref.  43,  x A2  =  0.5  x  (nocr3(^)s((y)s)5.  The  DSMC 
results  are  in  reasonably  good  agreements  with  the  results  of  Ref.  43,  although  a  noticeable  divergence  from 
correlation  values  occurs  at  low  values  of  the  scaling  parameter,  which  correspond  to  low  number  densities 
or  high  temperatures.  Note  that  the  empirical  correlation  underpredicted  the  data  in  that  region. 

V.  Water  cluster  size  distribution  in  nozzle  flows 

The  last  part  of  the  validation  and  numerical  analysis  of  the  presented  condensation  model  is  focused  on 
the  nucleation  and  evolution  of  small  water  clusters  in  a  small  conical  nozzle  nozzle.  The  study  was  prompted 
by  the  availability  of  quality  experimental  data44  on  terminal  size  distribution  of  water  clusters  in  the  range 
of  flow  conditions  where  the  Knudsen  number  was  relatively  large  so  that  the  computational  cost  of  using 
the  DSMC  method  is  not  prohibitive  (although  still  rather  high).  The  results  are  presented  for  the  nozzle 
geometry  of  Ref.  44.  The  nozzle  is  a  conical  nozzle  with  a  41°  opening  angle,  a  total  length  of  2  mm,  and  a 
throat  diameter  of  50  /im.  The  smallest  stagnation  pressure  considered  in  Ref.  44,  1.577  bar,  is  used  in  this 
work,  with  the  corresponding  stagnation  temperature  of  495  K.  Two  constant  temperatures  of  the  nozzle 
surface  were  considered  in  the  present  computations  to  study  the  effect  of  surface  temperature  distribution, 
the  baseline  495  K  (stagnation  temperature)  and  300  K  (room  temperature).  Since  the  background  pressure 
effect  in  the  experiment  is  believed  to  be  small,58  the  flow  expansion  into  the  vacuum  is  modeled. 

The  axisymmetric  capability  of  SMILE  was  used,  with  the  total  number  of  simulated  molecules  and 
collision  cells  about  80  million  and  10  million,  respectively.  A  uniform  400x100  grid  was  used  for  sampling 
of  macroparameters  and  distribution  functions.  The  Larsen-Borgnakke  model  with  temperature-dependent 
relaxation  numbers  was  used  for  energy  transfer  in  monomer-monomer  collisions,  and  the  reflection  of  par¬ 
ticles  on  the  nozzle  surface  was  assumed  to  be  fully  diffuse.  Uniform  inflow  conditions  were  imposed  at  the 
nozzle  throat,  calculated  from  the  isentropic  flow  relations  at  a  constant  specific  heat  ratio  of  1.4. 

The  first  set  of  results  presented  here  shows  the  effect  of  the  condensation  on  the  gas  flow  inside  the 
nozzle  and  in  the  plume  near  field.  The  gas  translational  temperature  and  axial  velocity  is  shown  in  Fig.  5 
for  two  cases,  the  baseline  condensation  model  and  the  condensation  turned  off.  The  results  show  that  the 
condensation  practically  does  not  change  the  flow  parameters  inside  the  boundary  layer.  This  is  expected, 
since  the  temperature  in  the  boundary  layer  is  higher  than  in  the  coreflow,  and  the  nucleation  near  the 
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surface  are  not  likely.  Near  the  centerline,  though,  the  condensation,  being  an  exothermic  process,  results 
in  a  significant  heat  release.  The  temperature  in  that  region  visibly  increases,  with  the  maximum  change  of 
over  20  K.  The  higher  temperature  for  the  condensing  flow,  accompanied  by  fast  translational  relaxation, 
causes  an  increase  in  the  axial  velocity  in  the  coreflow  of  about  30  m/s.  The  velocity  change  inside  the 
boundary  layer  is  negligibly  small. 


Figure  5.  Impact  of  the  condensation  on  the  water  translational  temperature  (K),  left,  and  axial  velocity 
(m/s),  right. 


Consider  now  the  effect  of  the  nozzle  surface  temperature  in  the  condensing  water  flow,  illustrated  in 
Fig.  6.  Generally,  a  higher  wall  temperature  should  result  in  thicker  boundary  layer,  which  is  consistent 
with  the  temperature  distributions  shown  in  Fig.  6  (left).  The  influence  of  the  wall  temperature  on  the 
coreflow  is  negligible  in  the  first  quarter  of  the  nozzle,  and  still  relatively  small  further  downstream.  The 
temperature  in  the  plume  coreflow  is  only  about  7  K  higher  for  495  K  case.  It  is  interesting  that  there  is 
practically  no  impact  of  the  wall  temperature  on  the  condensation  process  in  the  nozzle  core  flow,  as  the 
water  cluster  mole  fraction  given  in  Fig.  6  (right)  illustrates.  Although  the  region  where  the  condensation 
occurs  is  larger  for  the  colder  wall  case,  the  mole  fraction  profiles  along  the  nozzle  axis  for  300  K  and  495  K 
are  within  the  statistical  accuracy  of  the  calculations.  Generally,  the  cluster  concentration  increases  rapidly 
in  the  first  100  /irn  from  the  nozzle  throat,  and  reaches  its  maximum  of  about  1.5%  at  ~  200/im.  Then,  the 
evaporation  and  cluster  coalescence  result  in  some  decrease  of  the  cluster  concentration. 
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Figure  6.  Impact  of  the  nozzle  surface  temperature  on  the  water  translational  temperature  (K),  left,  and 
cluster  mole  fraction,  right. 


The  computations  conducted  for  a  higher  stagnation  pressure  of  5.177  bar  have  shown  relatively  little 
effect  of  pressure  on  both  gas  temperatures  and  number  of  clusters  in  the  flow  (see  Fig.  7).  The  boundary 
layer  thickness  decreases  with  increasing  pressure,  whereas  the  temperatures  in  the  coreflow  are  practically 
the  same  in  the  first  quarter  of  the  nozzle.  Further  downstream,  in  part  du  to  the  condensation  heat  release, 
the  temperature  in  the  coreflow  is  a  few  degrees  higher  for  the  larger  pressure  case.  The  cluster  mole  fraction 
is  somewhat  higher  for  2.177  bar,  but  generally  the  increase  is  less  than  ten  percent.  Again,  there  is  a  clear 
maximum  in  the  cluster  mole  fraction,  and  its  location  does  not  change  with  pressure. 

In  addition  to  a  small  increase  in  the  number  of  clusters  with  stagnation  pressure,  the  average  cluster 
size  also  increases  with  higher  pressure.  The  average  cluster  size  as  a  function  of  the  axial  distance  along 
the  nozzle  axis,  calculated  as  the  mean  value  of  monomers  in  a  cluster,  is  presented  in  Fig.  8.  The  number 
of  clusters  larger  than  dimers  is  larger  for  2.177  bar  case  at  any  point  in  the  flow,  and  the  relative  difference 
increases  with  the  distance  from  the  throat.  Due  to  larger  number  of  collisions,  it  takes  longer  for  the  higher 
pressure  case  to  reach  the  terminal  size  distribution,  that  appears  to  be  further  downstream  that  the  4  mm 
distance  from  the  throat  computed  in  this  work. 
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Figure  7.  Impact  of  the  stagnation  pressure  on  the  water  translational  temperature  (K),  left,  and  cluster  mole 
fraction,  right. 


Figure  8.  Average  cluster  size  along  the  nozzle  axis  for  two  different  pressures. 


Following  the  analysis  of  the  impact  of  different  physical  flow  parameters  presented  above,  a  sensitivity 
study  has  been  conducted  to  estimate  possible  effect  of  parameters  of  the  model.  In  addition  to  the  baseline 
computation,  the  cases  with  a  larger  Z  =  1  (compared  to  the  baseline  of  Z  =  6)  and  smaller  values  of  the 
cluster  evaporation  energies  and  cluster  heat  capacities  (both  decreased  by  20%)  were  used.  Note  that  these 
are  relatively  large  changes  that  are  believed  to  represent  the  limiting  cases  for  the  corresponding  parameters. 
The  cluster  mole  fractions  for  these  four  cases  are  presented  in  Fig.  9.  Note  that  the  clusters  represent  a 
trace  species  with  typical  mole  fractions  on  the  order  of  one  percent,  and  statistical  scatter  is  visible  for  them 
(the  species  weights  were  not  used  here  in  order  to  avoid  any  potential  energy  non-conservation  issues  that 
they  may  cause).  The  cluster  mole  fraction  profiles  were  therefore  smoothed  using  a  B-spline  approximation 
routine  in  SMILE.  It  is  interesting  to  note  that  the  mole  fraction  terminal  value  is  only  weakly  sensitive  to 
the  varied  parameters.  It  is  about  1%  for  all  cases  under  consideration.  The  transient  behavior  in  the  firs 
half  of  the  nozzle  significantly  differs  for  these  cases,  though.  The  decrease  in  the  heat  capacity  is  in  effect 
the  decrease  in  the  number  of  cluster  internal  degrees  of  freedom.  Smaller  number  of  internal  degrees  of 
freedom  of  a  cluster  means  that  when  the  Larsen-Borgnakke  routine  is  used  for  energy  redistribution  (the 
energy  is  split  proportionally  to  the  number  of  degrees  of  freedom),  it  will  take  more  collisions  to  raise  the 
cluster’s  internal  energy  to  the  values  exceeding  the  evaporation  threshold.  As  a  result,  more  clusters  are 
present  in  the  flow.  When  the  cluster  evaporation  energy  values  are  reduced,  smaller  number  of  clusters  are 
observed  in  the  flow,  since  they  become  easier  to  evaporate.  The  cluster  internal  energy  relaxation  number 
has  relatively  small  effect  on  cluster  populations. 

The  above  arguments  for  the  impact  of  different  parameters  of  the  model  are  still  applicable  when  the 
average  cluster  size  is  examined.  Only  the  change  in  Z  has  relatively  minor  effect  in  this  case.  To  the 
contrary,  the  decrease  in  the  heat  capacity  significantly  increases,  and  in  evaporation  energy,  increases,  the 
average  cluster  size.  Note  also  that  changing  all  three  parameters  seems  to  increase  the  relaxation  length  as 
compared  to  the  baseline  case,  where  the  cluster  fraction  and  size  distribution  come  to  their  terminal  values 
soon  after  the  nozzle  exit. 
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Figure  9.  Cluster  mole  fraction,  right,  and  average  cluster  size,  left,  along  the  nozzle  axis  for  different 
parameters  of  the  condensation  model. 


Comparison  of  the  terminal  size  distribution  obtained  for  the  baseline  model,  with  experimental  re- 
suits44,59  is  shown  in  Fig.  10.  In  the  experiments,  the  dimers  were  below  the  detection  threshold,  the  trimer 
population  may  be  somewhat  affected  by  that  threshold,  and  all  larger  clusters  are  believed  to  be  recorded 
without  significant  distortions.  For  side-by-side  comparison,  the  experimental  size  distribution  was  normal¬ 
ized  so  that  it  has  the  same  fraction  of  clusters  larger  than  dimers  as  in  the  computation.  The  results  show 
that  while  there  is  a  reasonable  agreement  between  the  numerical  modeling  and  measurements,  there  is  a 
maximum  for  the  cluster  sizes  of  5  and  6  observed  in  the  computations,  that  were  not  recorded  experimen¬ 
tally.  The  reason  for  this  maximum  is  the  corresponding  maximum  in  the  evaporation  energies  for  these 
cluster  sizes  (see  Table  G).  It  is  possible  that  more  accurate  or  general  values  of  evaporation  energies  (only 
one  water  molecule  orientation  was  considered  in  Ref.  50)  have  to  be  used  to  obtain  better  agreement  with 
the  data. 


Figure  10. 


Comparison  of  the  computed  size  distribution  with  experimental  data  of  Ref.  44. 


VI.  Conclusions 

A  first-principle  condensation  model  applicable  for  the  DSMC  method  is  constructed.  It  was  validated 
through  comparison  with  available  theoretical  and  experimental  data  on  condensation  rates  in  a  thermal 
bath,  dimer  mole  fractions  in  orifice  expansions,  and  cluster  size  distributions  in  nozzle  flows.  The  model  is 
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based  on  a  DSMC  model  of  recombination  reaction  for  the  collision  based  dimer  formation,  and  the  RRK 
model  for  the  cluster  evaporation.  Cluster  growth  is  modeled  through  cluster-monomer  and  cluster-cluster 
collisions.  The  energy  transfer  in  these  collisions  is  calculated  using  the  extended  Larsen-Borgnakke  principle. 
An  important  parameter  that  controls  the  rate  of  internal  energy  transfer  in  cluster-monomer  collisions  is 
analyzed,  and  the  sensitivity  of  evaporation  rates  to  this  parameter  is  shown. 

Two  gases  are  considered  in  this  work,  argon  and  water.  For  the  thermal  bath  relaxation,  the  present 
model  was  found  to  capture  the  equilibrium  constants  for  water  and  argon  and  nucleation  rates  for  argon 
fairly  well  in  the  considered  range  of  temperatures  from  100  K  to  350  K.  The  computed  slope  was  somewhat 
smaller  than  theoretical  for  lower  temperatures,  which  may  be  attributed  to  the  accuracy  of  used  physical 
and  numerical  parameters  of  the  condensation  model.  Comparison  of  terminal  mole  fractions  in  a  sonic 
orifice  expansion  dimer  with  the  available  semi-empirical  correlations  showed  that  the  new  model  agrees 
reasonably  with  the  correlations,  overpredicting  the  latter  by  about  a  factor  of  two  for  lower  pressures  and 
higher  temperatures. 

Sensitivity  of  the  clusterization  to  the  cluster  internal  energy  relaxation  number,  cluster  evaporation 
energies,  and  heat  capacities,  has  been  examined  and  shown  for  water  expansion  through  a  conical  nozzle. 
Additional  studies  are  necessary  to  further  refine  the  utilized  values  of  the  above  parameters.  Comparison 
of  the  cluster  size  distribution  with  the  available  experimental  data  has  been  conducted,  and  reasonable 
agreement  between  the  computations  and  the  experiments  have  been  shown. 
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